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Abstract 

Single-layer carbon, or graphene, demonstrates amazing transport properties, such as the min- 
imum conductivity near ^ independent of shapes and mobility of samples. This indicates there 
exist some unusual effects due to specific Dirac dispersion relation of fermion in two dimensions. 
By deriving fermion-lattice interaction Hamiltonian we show that Berry phases can be produced in 
fermion states around two Dirac points by relative rotations of two sublattices. The Berry phases 
in turn remove the degeneracies of energies for states near the Fermi surface, leading to a dynamical 
instability of the lattice with respect to the rotations. By considering the Berry phases emerging 
in an uncertain way on fermion wavefunctions in vicinities of the Fermi surface, the conductivity is 
calculated by using the Landauer-Biitticker formula together with the transfer-matrix technique, 
verifying ~ quantized minimum conductivity as observed in experiments independent of shapes 
and sizes. The relationship between the chaotic structure of fermions due to the Berry phases and 
the classical transport properties are discussed. The physical meaning is profound as this relation- 
ship provides an excellent example to elucidate the mechanism of quantum-classical transition. 
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As a potential candidate material for carbon-based electronics, the graphene, a two- 
dimensional (2D) honeycomb lattice of carbon atoms, has attracted much attention 
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dispersion relation of the massless Dirac fermions in graphene [3, [J The Fermi level of the 
system can be varied by changing the gate voltage. Particularly, there exists a universal 
maximal resistivity of the order of 6.5kf2 for all samples with the Fermi level near the 
Dirac points, independent of their mobility. By increasing the Fermi level the conductivity 
linearly increases j^]. The value of the minimal resistivity {-P^ ~ 6.5kfi) [3] implies the 
existence of at least 4 channels (including spin degrees of freedom) in a square of unit area 
and the number of channels increases with increasing the sample width. Although there 
are 4 ballistic channels at a node of the Dirac dispersion relation in a graphene lattice, 
the experimental fact of constant conductivity is inconsistent with the quantum-mechanical 
description from which all states are localized by extremely weak disorder in 2D [9]. The 
non-zero minimal conductivity has been investigated by several theoretical studies using 



different methods 



[lq, 



11 



12 



13] . The minimum conductivity was correctly demonstrated, 
but as to the reason for it and as to its exact value being e 2 / h or e 2 / nh there is no consensus 
at the present. By carefully examining the interaction between Dirac fermions and lattice 
motions, we show that the Berry phases produced by relative rotations of two sublattices 
can induce dynamical instability and cause constant conductivity 

The main features of the band structure in graphene can be well described by a tight- 
binding Hamiltonian with one tt orbital per site on the graphene lattice Q: 

H = W(<4 a n' + a^,a n ), (1) 

(n,n') 

where cv n (a n ) creates (annihilates) an electron at site n, (. . .) denotes the nearest-neighbor 
(NN) sites, and t nn > is the NN hopping. Here, we omit the diagonal terms by shifting the 
energy zero as no diagonal disorder is considered, and the spin indices are not explicitly 
included. The hopping integral t nn > depends on distance between NN sites. Considering a 
relative displacement d = (d x , d y ) between two sublattices as shown in Fig. 1, the lengths 
of three bonds in a cell become 



k = J[dsm{6 - a,)] 2 + [l - dcos{6 - c^)] 2 , (2) 
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where 9 is the angle of d related to x axis, 1$ is the original bond length, and a>i = 0, — ^- 
for i = 1,2,3. For small displacement, we only keep the first-order terms of d and obtain 
k ~ — dcos(9 — otj). Then the corresponding hopping integrals can be written as 

£j = *o + Xdcos(9 — oij), (3) 

where to is the original hopping and A is a coefficient describing the linear dependence of U 
on li in the case of small d. On the other hand, the elastic energy trends to keep d = 0. The 
in-plane displacements occur in two vibrational modes along x and y directions. We adopt 
variables d x and d y to express these two vibration modes. Using the Bloch transformation 
for graphene lattice and keeping terms of the first orders of d x , d y , and electron wave vector 
k related to the Dirac points, the Hamiltonian becomes 

H = - [(toky? z + Xd x l) ®a x - {t k x l - \d y T z ) <g> a y , (4) 
1 k 

where 1 and f z are unit and Pauli matrices acting on two irreducible Dirac points, and a XjV 
are Pauli matrices on two sublattices. So the dispersion relations around two Dirac points 
are 

3 3 

= i- \t Q {ky + ik X ) + \{d X ~ idy)\ , = i- \t Q (-ky + ik X ) + \{d X + idy)\ . (5) 



It is interesting to note that the diabolical points of the dispersion relations are moved by 
sublattice displacements. The dependence of energies on real d for a given k also has 
a diabolical point determined by (d x = —t k y /X,d y = tok x /X) for E^ and by (d x = 
toky/X,d y = —t k x /X) for E^ as shown in Fig. 2. 

From the Berry theorem a circle of rotation in a parameter space can cause a Berry phase 
±7r in a wavefunction whose eigenenergy has a diabolical point enclosed in this circle [bj) ]. 
If there exists a circular vibration with amplitude df = t kf/X with kf being the magnitude 
of the Fermi wave vector, in every circle of the vibration the states inside the Fermi circle 
change sign due to the acquired Berry phases but the states outside the Fermi circle do not. 
Now we consider an inside state j^i) and an outside state \ip2) situated in an infinitesimal 
vicinity of a k point at the Fermi circle. They can be degenerate because of the continuity 
of the band. So we can construct a 4>+-ip- basis to represent these two states: 

llM = ^(lV>i) ± l^2». (6) 
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FIG. 1: Sublattices A (red points) and B (blue points) on graphene lattice. With relative displace- 
ment d, sublattice A is moved to positions marked with light magenta circles. 

The asymmetric acquirement of a ±7r Berry phase in IV'i^) during a circle of circular vibra- 
tion df leads to a hybridization of ^>±(k) and removes the degeneracy. This can be described 
by an effective Hamiltonian 



where 6f is the Fermi energy and g is the hybridization strength depending on the rotation 
rates. Hamiltonian Hb raises one level to e/ + g and lowers the other to e/ — g. Since 
the former is empty and the latter is occupied by fermion, this leads to the lowering of 
the fermion energy by g. Thus, the energy gain per cell of the fermion system due to the 
vibration df is E g ~ li^k^gD where D is the width of a narrow stripe around the Fermi 
circle in which the fermion energies are affected by g. At the same time the energy per cell 
of the circular vibration df is proportional to d/j and can be written as E v ~ cd\ with c 
being a prefactor. So, if E v < E g , viz. 



H B =Y, \^)e f (^a\ + g(\i/>+)(i/>-\ + H.c), 



(7) 



a=± 



kf < k v = 



2-KgDX 2 
ct 2 



(8) 
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FIG. 2: Ei as functions of d for different values of k. 

a dynamical instability occurs with undamped fluctuations of rotations and Berry phases. 

The condition (jSJ) for the dynamical instability can be easily satisfied since k v is finite 
and positive while kf — > for the Fermi level situated at the Dirac points. It could still be 
satisfied even by slightly shifting the Fermi level with the gate voltage V g . The value of k v 
can be estimated as follows: From the Harrison formula, the hopping integral of NN sites n 
and n is t nn > = t ( j- 2 - 



■' rt n 



15] . So A ~ j 2 - ~ 1.4to A -1 . c is related to the elastic coefficient 
s e of the C-C bonds by c ~ |s e , and s e ~ 43.37eV A -2 [16]. (7 is related to the vibration 
frequency ujq by g ~ frc^o, and D is mainly due to the uncertainty of the Fermi energy in 
the vibration and can be expressed as W ~ ^f^- From the value of s e and the mass of 
carbon atom we have Tiujq ~ O.leV. All these can give k v ~ 10 -3 in units of Iq . This value 
is considerable as the whole Dirac dispersion relation is valid for small k. 

The circular vibrations could also propagate in the space with different wavelengths longer 
than the lattice spacing (long-wavelength limit) and cause mesoscopic inhomogeneity of the 
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created Berry phases. The spatial and temporal uncertainty of Berry phases attached on 
fermion states by such circular vibrations has radical effects on the basic structure of fermions 
in graphene. To illustrate this we adopt a model with random signs of NN hoppings t nn i = t 
or — 1 to mimic the fluctuations of Berry phases ±7r and calculate the conductivity for a 
rectangular sheet with different width W and length L connecting to two reservoirs with 
Fermi energy at the Dirac points (Ef = 0) by using the Landauer-Biittiker formula: 

2£e 2 v 1 

° Wh Vcosh 2 ( 7z L)' { > 

where 7/ is the Lyapunov exponent of the Zth channel in the system with width W calculated 
with the transfer-matrix technique. The prefactor 2 is from contributions of two spins. The 
obtained conductivity as a function of W/L for different sizes is shown in Fig. 3. Except 
for small values of W/L, the conductivity has almost a constant value near independent 
of W/L and W, in consistence with the experimental findings. Surprisingly, this value is 
even larger than ^ obtained in ballistic graphene 17 1 in spite of the disorder introduced 
by random signs of hoppings. In a perfect 2D Dirac fermion gas, the difference between two 
successive Lyapunov exponents in Eq. is £7 = 2k/ W due to the Dirac fermion dispersion 
relation [171 ] . In the present model the values of W5j are distributed within [1, 2.5], smaller 
than that in the perfect case by a factor about 7T, leading to the 7r times of the conductivity 
as can be derived from the expansion of Eq. (191) . Thus, the 7r factor difference in the 
conductivity is not trivial and reflects an essential change of the fermion properties due to 
the Berry-phase induced instability. 

Besides the phase fluctuation, there are complicated level crossings in the dependence of 
energy on k and d. The energy degeneracies at crossings can be easily modified by slight 
disturbing, resulting in chaotic energy structure. Thus, the global quantum-mechanical 
description of the Dirac fermions could be modified with a semi-classical scenario: The 
whole sheet can be divided into areas of size &o which are classically connected together 
as there is no quantum coherence among them. This is equivalent to an ohmic network 
of small conductors each of which represents an area of size bo. Thus, the conductivity of 
the whole system is determined by the conductivity of every area which can be quantized. 
For the Fermi energy Ef = there are only two transmission channels in an area providing 
minimum conductivity a — By applying the gate voltage V g the Fermi energy is shifted 
as Ef = f3V g with j3 being a prefactor and the number of transmission channels in an area is 



2/3 3Kto° ' P r °ducing a conductivity a ~ • The experiments show a linear dependence 

of a on the gate voltage V g with a finite slope, suggesting a finite coherent size b . From the 
slope one could estimate ~ 0.2y _1 . 

In summary, we show the specific coupling between the Dirac fermions and the circular 
vibrations in graphene which can lead to a dynamical instability due to the produced Berry 
phases. The calculations from this scenario provide a constant conductivity consistent with 
the experiments. The results indicate that a picture of an ohmic network consisting of small 
quantum conductors is suitable to account for the peculiar transport properties in graphene. 
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FIG. 3: Calculated conductivity as a function of ratio W/L for rectangular sheets with different 
transverse sizes W. Inset in the lower panel: Distribution function P(W5'y) of W5j. The green 
dashed line indicates the 5-function distribution at 2tt in the case of perfect sheet. Units of W and 
L are cell sizes of graphene lattice. 
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